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ABSTRACT 


A discussion of s2veral non-numerical techniques that 
are useful in FORTRAN programming is presented. The use of 
these techaigues is than illustrated with a problen oriented 
language called CAL-NPS. This last program is a derivation 
9f a code named CAL written by Professor EB. L. Wilson of the 
Gniversity of California, Berkel2y, California. 
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I. INTRODUCTION 


In writing Problem Jriented Computer Languages several 
non numerical applications of computer programming must be 
mastered and the literature about such activities is nearly 
nil for the PORTRAN language. It is to these problems that 
this work addresses itself. The professional programmer 
would probably elect to write these codes in machine 
language but the user who only knows FORTRAN cannot afford 
tine to learn the various machine languages in use today. 
All of the features discussed pb2low are progranmable in 
FORTRAN ani illustrations are given in each use. 

The CAL languag2? is then -‘us2di to illustrate how these 
features are incorporated in a problem oriented language. 
[The CAL program was jisveloped by Professor Edward L. Wilson 
of the University of Salifornia, Berkeley, Caltix.rnia, LOL 
Structural Mechanics and Structural Engineering. 

CAL conubines matcix manipulation routines with direct 
stiffness computation options to produce a program for the 
automated analysis of structures. CAL also has significant 
Capapiliti2s as an instructional td0l in linear algebra. 

In orier to use any structural analysis program it is 
necessary to idealize the structure into a series of joints 
connected by structural members. The joints are commonly 
referred t> as nodal points andthe structural members as 
elements. The program input consists of: 

a) nodal point locations; 

b) information on how nodes are connected; 

>) boundary conditions that are applied to the 
nodes; 

d) properties of the elsnents; 

2) the loading to be analyzed. 


The progran output, at a minimnun, consists of nodal point 





deflectioa values and normally includes nodal forces and/or 
element stc2ss valués. CAL breaks this procedure down into 
a series of simple steps under th2 control of the user. It 
was designed to execute rapidly on small problems, so 
Students can quickly see results of an analysis. There are 
options which allow the user to d2bug data without printing 
praviously obtained results. Looping operations are also 
available. This allows a user to program iterative 
numerical algorithms, greatly expanding the usefulness of 
the program. In this version, CAL will solve problems with 
about 50 degrees of fraeiom (DOF). 

The orimary aim of the thesis work presented here was 
to modify CAL for use on the IBM 360/67 computer at the 
Naval Postjraduate Schsol (NPS). The NPS version of CAL has 
both interactive and batch operating modes, and is nearly 
machine independent. | 

Most of this thesis effort was not involved in the 
numerical procedures ased for calsulations, but rather in 
the non-nuaecical areas of data input and processing. fFor 
this reason, numerical procedures will not be discussed. 
Subsequent chapters will provide: 

a4) a substantive discussion of some non-numerical 
technigues useful in writing scientific programs such as 
SAL ; 

b) a discussion of the internal organization of 
~AL-NPS; 

Sc) instructions for the implementation of CAL on a 
different computer. 
The appendices contain a listing of the FORTRAN code and an 
Instruction Manual for CAL-NPS. 

Modifications to the progras were carfully made to 
avoid changes in the instruction nanual written by Professor 


Wilson. f2n commands have been adied and three commands 





have been 2xtended for more flexibility in the interactive 
use of th2 program. The original instruction manual has 
bean modified to reflect these changes and is reproduced in 
Appendix 3. The author gratefully acknowledges Professor 
Hilson's parmission to use the manual. 

The noo-~numerical techniques discussed in Chapter II 
are this author's effort to explain how several 
state-of-the-art programming techniques work. It is hoped 
that readers will find this section beneficial in their own 
programming efforts. Th2 CAL organization discussed in 
Chapter III is an excallent outline for general scientific 
programs. Additionaly, Chapter III provides details for 
writing user supplied subroutines. For readers who are 
interestel only in how to use CAL, Chapter II and IV may be 
Di Gee d : 





If. NON-NUMERICAL TECHNIQUES IN FORTRAN 


Chapter II 15 used to discuss FORTRAN programming 
techniques for writinj general scientific programs. The 
reader is assumed to have some faniliarity with the FORTRAN 
language and computer programming (1.e., completed cCS2700 
FORTRAN Programming or equivalent at NPS). fter 
loatroductorcy remarks on the attributes of a scientific 
prodgran, idiynamic dimensioning is discussed in Section A, 
fr2e and odbject time formats are discussed in Section B, 
data management is discussed in Section C, and finally, 
Paepang 25 discussed in Section OD. The terms will be 
defined as encountered in these s@stions. 

A gen2ral purpose program for the computer solution of 
4 class 9£ scientific problems ‘should be simple to use, 
flexible ani reliable. This reguices a program that: 

a) 1s modular, i.e., sections of program may be 
deleted, inserted or modified with ease; 

b) has num2rical methods that are stable ani accurate 
for the tyoe of problem encountered; 

c) may reguire the user to have a knowledge of 
munecical methods necessary for the selection of the 
appropriate solution technigue or boundary conditions but 
do2s not reguire knowledge of computer programming beyond 
the ability to prepare the input data for a progran, 
initiate 2xecution of a program, and create data files; 

d) is not limited by array sizes in dimension or common 
statenents; 

e) permits iteration, i.e., the repeated execution of a 
user selected sequence of operations; 

f) ceyuires a minimum of c2antral processor time and 


internal computer memory space (core) for execution; 





A very useful method for meeting these requirements is 
to structare the program arouni the use of a  problen 
oriented language. Such a language consists of a set of 
mnemonics cepresentinjgy each possible step inthe general 
solution of the problen. 

A mnemonic is a word that serves as a memory aid to 
identify some thing or some action. In FORTRAN, the 
mnemonic 59S is used to represent the process of computing 
the cosine of an angl@. All mnemonics can be thought of as 
representing a function subprogram, which can be very 
complex. fhe use of Mnemonics to represent various 
Sp2rations in the formulation and solution of a particular 
problem is a very attractive methosi of writing a program for 
the general solution of a class of scientific problems. 
This technigue promotes modularity and greatly simplifies 
use of the program. The mnemonis used for a particular 
operation or computation should be representative of the 
Operation. For instance, a programmer might select the 
mnemonic SOLVE for th2 operation of computing the vector x 
from the aatrix equation Ax = B where A and B are known. 
[Th2 mnemdaic identifying a particular operation together 
With a list of arguments necessary for the computation is 
read from one data cari. This cari is followed by a series 
9f data cards to provide data, not currently in core 
Seotrage, £f9r computation. MThe ability to enter the 
information in free format or t9 allow object time format 
specification dramatically reduces the work required to use 
the progrin. 

References (1), (2), and (3) provide examples of 
problem orianted language programaing. Apart from these few 
examples, there is very little written on the non-numerical 
techniques needed to write a program for the general 


solution of a scientific problem using a problem oriented 


elO 





language for the solution. References (4) and (5) are 
useful 2xceptions and provide sone very valuable 
information. Th2 numerical techniques for execution time 
arcay dimensioning, execution time and free format, data 
management and the automated repetition of a user selected 
sequence of instructions will now be discussed in some 
detail. 


aA. DYNAMIC DIMENSIONING 


In classical programming, array storage is reserved 
either in common or in dimension statements. Either way, 
the dimensions are usually set to the largest values 
2nvisioned for the type of problem under solution. When 
Small or nedium sized problems are solved, a considerable 
amount of core space is reserved but not used. Avoiding the 
wasted space normally requires re-doing the dimension and/or 
common statements in the entire program and then 
re-compiling the whole progran. 

Compilation is the first step in preparing a program 
for execution. In the compile step, the FORTRAN statements 
are interpretei to produce a set of machine language 
instructions, called an object deck, that the computer is 
capable of following. The next step is to link-edit the 
object deck. The link-editor collects all of the machine 
language instructions needed to read data, write, and 
evaluate functions called by th2 program; puts then in 
logical order with the object deck; and prepares an index to 
locate various sections of the program. The output from the 
link-editor is called the load asddule. The final step is 
the actual execution of the progran. This is called the go 
step and consists of storing tha load module in core then 


executing the instructions in sequence. Pigure 1 


vt 
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demonstra:es this process graphically. Usually it is 
possible t> use a compiler that will optimize the machine 
language coie to winimize the run time for a program in 
execution. Optimizing compilers take more time and require 
more core space to compile than do non-optimizing compilers. 
It is possible to compile and link-edit a program then store 
the load madule ona mass storage device (magnetic tape or 
disk). This always results in saving central processor 
(CPU) tin2 for subsegient program runs and can result ina 
Significant reduction of the amount of core required. 


COMPILER 
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CARDS 


MODULE 





PROGRAM 
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Program Processing 


PIGURE 1 
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The aiternative to the classical style of programming 
1s dynamics dimensioning. [In ordar to understand the concept 
Sf dynamic dimensioning, itis necessary to have some 
aporeciationo for how a computer handles a common statement. 
Inside the computer, the common statement will cause the 
reservation of sufficient bytes of memory to store listed 
variables. A byte is a segment of a computer word which can 
contain either one alpha (hollerith) character or about two 
decimal jigits. If the same coanon requires a different 
number of bytes in different subroutines, the largest 
reguirement will be th2 one reserved. Any type of variable 
may occupy the reserv2d storage as long as the address of 
th2 starting byte of the variable is consistent with the 
type of variable. Different types.of variables frequently 
have a different number of bytes in a word. For instance, a 
complex variable will normally hav2 twice as many bytes per 
word as a ceal variable. To illustrate this point, consider 
a program where all variable names beginning with A are 
double precision (8 bytes), all variable names beginning 
with Z are complex double precision (16 bytes), and 
variables beginning with N are integers (4 bytes). The named 
common EXAMPLE might be listed as follows in different 


subroutin2s: 


COMMON /EXAMPLE/ N(12) (1) 
COMMON /EXAMPLE/ A (3,2) (2) 
COMMON /EXAMPLE/ Z(3) (3) 


In each case the same number of bytes of storage are 
reserved (48) but the common has a Significantly different 
type of variable! Introductory FORTRAN courses normally 
teach the student that a particular common should always be 


an exact duplicate each time it appears inthe progran 
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Without aach other discussion of the uses of common 
Statements. The example lines above contrast sharply with 
this cautious style but can prove very useful in actual 
problem solving. Table 1 gives th2 address of the starting 
byte for 2ach variable in EXAMPLE. 


BYTS VARIABLES 

1 N (1 Als Vy eect) 
ee 
E es. xs 
> pees Cay vie 20) 
25 N(7),A(1,2) 

$3 gat a> Z(3 
a7 Neate ¢ cunge v=) 
re N 1} (3,2 

45 N (12) 


Byte Aidresses for Variables 
ASCE als 


Meee that N(3), A(2,1) and the imaginary portion of 2Z(1) all 
have the same byte as their starting address. Additionally, 
48 bytes would have been reserved even if the array 
dimension of both (2) and (3) had been set at one since (1) 
Still requires 48 bytes. One more feature of importance is 
that A(2,1), A(1,2) ani A(3,2) contain the same information 
as the coresponding imaginary portions of the Z vector. This 
fact is guite useful in programming with complex variables. 
The common EXAMPLE above denonstrates how a computer 
stores a two dimensional array but amplification may prove 
helpful. The right most subscript in a multiply dimensioned 
arcay is always the one that increases Slowest as the array 
is going into comput2r storage. The address in memory of 
an element of a linear array 15 somputed much faster than 


the address of an element from an array with multiple 


ee ee ee ee ee eee ee 





SupScripts. The address of the element A(I,J) of a two 
dimensional array can be computed in a corresponding linear 
array as follows: 
L = ((J-1) *NR + I - 1)*IPR + 1 (4) 
wher=2 

NR is the number of rows in the matrix 

IPR is the ratio of bytes per word of the matrix to 
bytes per word of the linear array. 
For example, compute the index in WN of A(3,2) from the 
common EXAMPLE: 

here I=3, NR=3, J=2 and IPR = 8/4 = 2 

L = ((2-1)*3 + 3 -‘1)*2 + 1 = 11 
Sompare this answer with the results in table 1. 

Dynanic dimensioning is a process that utilizes the 
computer's ability to store different types of variables in 
the same common array. Array storage is reserved in blank 
common with a larg2, one-dimensional, integer array 
(hereafter called th2 main array) for all subscripted 
variables needed for a particular computation step. This is 
normally ione in a very short main program. Problem size can 
b2 controlled by the dimension of the main array in the main 
program. The dimension of the main array in all subroutines 
is set az 1and never changed. Typically, to change a 
problem size, two cards of a six card main program need to 
be alterei. Seo figure 2 for an example. The main progran 
ne2ds to be compiled, but there is no need to use an 


Optimizing compiler. Th2 remainder of the program can be in 


2 





a load moijule on a mass storage device. 


MAIN PROG 
SIMMON MAX,L (6000) 
MAX = 6000 
CALL DRIVER 
STOP 
END 

C----- SUBPROGRAM TO CONTROL EXECUTION 
SUBROUTINE DRIVER 
=OMMON MAX,L(1) 
SALL COMPUTE(L(N1),L(N2),...,N3,N4) 
RETURN 
END 

C-----SUBROUTINE FOR COMPUTATION 
SUBROUTINE COMPUTE (A,B,-..,NROW,ICOL) 
DIMENSION A(NROW, 1) ,B(1) 
RETURN 
END 


Dynamic Dimensioning Sample Progran 


FIGORE 2 


The main program calls a subprogram which controls 
execution, reserves storage in the Main array for 
subscriptei variables, and provides the address of array 
arguments in the main array to the computation subroutines. 
This subprogram must also ensure that there is sufficient 
room inthe main array before a subscripted variable is 
allowed t> be stored. Boundary misalignment occurs if the 
type of vaciable being stored requires more than one integer 
word for storage, and the Starting address of the 
subscripted variable in the main array is not an address 
that is compatible with allowable addresses for the variable 


being stoced. Some compiliers automatically correct for 





boundary alignment but often this iS accomplished with a 
considerable p2analty in the time required for execution. 
Consider, for example, that a real double precision variable 
On an 134 360 computer requires two integer words for 
storage. A boundary misalignment will occur unless thea 
double precision variable is at an odd integer word address 
counting from the first integer word in blank common. In 
pguee 2, MAX 1S the first inteajges word and L(1) as the 
second word in blank common. L(2) is, therefore, the first 
address available for the storage of the double precision 
element. L(2) and L(3) are requir2i for the storage of one 
doable pr2tision number. If the double precision variable 
Beatted 12 h(i), Lt» would be out of proper alignment for 
this type of variable. A Simple way to compute a proper 
Starting address is: 

N1 = NSIZE + MOD( NVAR + NSIZE , IPR) + 1 

wher: 

N1 is the starting addr2ss to be computed. 

NVAR is the number of integer words used by the 
Vaciables in blank common preceeding the main array, L. 

NSIZE is th2 number of words in the main array 
that have been used. 

IPR is the number of integer words required to 
store the variable being processed. 

MwOoDeI,J) is ane NST Standard function which 
returns th2 remainder from (I/ J). 

A method to keep track of the starting addresses for 
Subscriptei variables will be discussed later. The amount 
of storage available may be compared with the amount 
required for storage of a subscripted variable as follows: 

LEFT = MAX ~ (N11 + NEL * IPP ) + 17 

where 


LEFT is the number of integer words remaining. 


Ly 











MAX is the number of integer words in the main 
arcay. 

N1 as the Starting address in the main array for 
the variable. 

NEL is the number of elements in the variable. 
Note that when blanks are left between the end of an array 
and the stact of the next array t> compensate for boundary 
alignment, NEL of the leading array is effectively changed. 

IPR 1s the number of integer words required per 
element. 

Posizive steps must be taken to prevent LEFT fron 
becoming negative. One method to compute the number of 
words that have been used in the main array is 

NSTZE = BAX = LEFT 

wher? all variables are used as defined above. 

Figuc2 2 is an exanple of some programming necessary to 
imolement the dynamic dimensioning concept. The variable 
MAX and the dimension of the L array must have the same 
value in the main program. Outside the main program, the L 
arcay ne24ds only to be dimensioned one since adequate 
storage is already reserved. In the control subprogran, N1 
and N2 are used to provide starting addresses for 
subscripted variables to the computation subroutine. N3 and 
N4 are usei to provide array dimensions. The dimension 
Statement in the compatation subroutines can use variables 
for array subscripts. Linear arrays and the final subscript 
of multiply subscripted arrays need only be set to one as 
shown in figure 2. 
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B. FREE AND OBJECT TIME FORMATS 


Recall that in classical programming, data entered into 
the computer must be entered exactly in accordance with the 
format statement that has been compiled with the prograa. 
For instaice, suppose that the coordinates and a reference 
nun ber are required for point in 3-s pace in 
mor mat (3F15.5,1510). Pigure 3 is the data card required for 
the point (1.0,-1.0,1.0) assigned reference number 3. Note 
that if ta2 3 is not in card column 55, the reference number 
assigned t> the point will be wroog. When the input device 
1s a terninal, the probability of making an error on 
entering data with this type of fixed forgat 1s ‘Hagh. 
Errors in entering data frequently require the user to 
restart th® program and re-enter all the data to correct a 


v 


Small mistake. 
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: 
Sample Fixed Foraat Data 


FIGURE 3 
Jbjest time format (also called variable format) 


permits th2 user to enter the format for subsequent input or 
output operations as a part of the data deck at the time the 
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program enters execution. This means that it is possible to 
read an array prepared in (F10.0) format and an array 
prepared in (5E15.7) format without changing the format 
‘statement and recompiling the progran. 

Object time formatting is accomplished by reading a 
card containing the format for the data to follow, then 
reading or writing the data as specified in the format. A 


Sample of the reguired programming to read data is: 


DIMENSION FOR(20),A(3), «~ 


100 READ (3,200) FOR 
200 PORMAT (20A4) 
300 READ(5,FOR) (A(I) ,1=1,3) ,0 


Statement 100 represents an ordinary read command in FORTRAN 
according to the format statement numbered 200. This format 
Statement causes 80 solumns of hollerith data to be read (4 
columns er1ch into 20 variables). Statement 300 causes 4 real 
rumbers to be read into the vector A. However, instead of 
being read accordiny to a nuabered format statement 
contained in the program, these numbers are read according 
to the information contained in thea variable FOR. FOR might 
Sontain: 
(3F4.0,1I1) 

CauSing the coordinates of the point discussed above to be 
read as the vector A in format (3F4%.0) and the reference 
number to be read as the integer J in format (I1). Pigure 4 


Shows the cards reguired to enter the point and reference 
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Sample dbject Time Dimension Data 


FIGURE 4 


The first non-blank character on the data card containing 
the format must be an open parenthesis. The last non-blank 
character gaust bea close parenthesis. The characters in 
between must constitute a valid format for the data to be 
read or written. 

Object time dimensioning allows the user to adjust 
formats to suit the data being input. However, if the user 
nejylects to enter a format or 2nters one containing an 
error, the resulting problems can be frustrating. Free 
format relieves the user of the difficulties associated with 
auSing any type of predefined format. Figure 5 is an example 
of the data card required to enter the point and reference 


number discussed earlier. This example uses a conma (,) as 
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an end of field symbol and a blank () as the end of record 
symbol. 
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FIGURE 5 


The first step of programming free-forBat input 2s to 
cread the data card into an array called the card-image 
arcay. The data card must be read in format (80A1) so that 
the Jthn element of the card-image array represents the 
character in column J of the data card. It is important to 
note that the value stored in the computer for a digit read 
as a hollerith charast2r is different from the value stored 
for the sane digit read as a number. The card-image array 
is then processed element by element to form the input 
variables. These input variables may be integers, real 
numbers oc character strings, however, the variables must 
still be in the order expected by the program. If the 
expected input is a data card containing a step name and 
matrix argiment list, a card containing real numbers will 
Cause an error. Free format input allows the program and 
not the system error handling routine to determine 


appropriate corrective action. 
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The first step in processing an element of the 
Card-imag> array is to identify the character storei in the 
el2ment. IdWw precise the identification must be depends on 
the application. The character being used as an end of 
field synbol and the character being used as an end of 
record mist always De identified precisely. Pep for 
instance, the expected data 1S a step name and matrix 
argument list, then the programner might chose to identify 
all alphabetic characters as one type of synbol, all digits 
aS a seconi type of symbol and all characters not otherwise 
identifiei asa thiri type of symbol. However, if the 
expected data contains a real humber that may or may not 
have an exponent the choice of characters nights be; "0.122, 
emer ye 06CGl UU OD, UhylU te eg D,lCUE,6CUCUtmrnd «€(Others. The Wariable 
NSYM will be used to indicate the type of character. 

Line Hastory of the processing is important in deciding 
what to 40 with a character. For instance, in the 
assemblinjy of the rceal number 1.21, there are two types of 
character, digits and a decimal point. If the character 
encountered is a digit, the proper assembly action depends 
90 whether the d2cimal point has occurred. If the decimal 
point has not occurred, the number being assembled must be 
multiplied by 10 and then the digit added to the number. On 
the other hand, if the decimal point has occurred, the digit 
must be divided by 19 then added to the number being 
assembled. The power of 10, N, is a counter representing 
the number of places to the right of the decimal this digit 
belongs. fhe varialbe IACT will be used to indicate the 
appropriate assembly action. 

In adiition to integral and fractional digits, a nuseber 
may carry an algebraic sign and/or a2 sign in the exponent. 
The sign aay occur at the beginning of the number or at the 


start of the exponent, however, a sign encountered elsewhere 
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Wall definetdy be an W@rror. It follows that for the 
assembly process of real numbers, the program performs one 
9f the following operations: (a) sign identification, (b) 
integral digit assembly, (c) deztinal point identification, 
(d) fractional digit assembly, and (e) exponent assembly. 
These prosesses are refered to as states. If the state is 
integral ligit assembly and the character encountered is a 
Sign, then obviously an error has occurred. On the other 
kand, if the character encount2r2ei is a decimal point, then 
there 1S 19 error but the state for the next operation needs 
ts be changed to fractional digit assembly. A variable named 
BanpaATe 1s introduced to indicate which one of the above 
operations the program last performed. Now, continuing with 
the assemoly of areal number, the variable ISTATE may be 


assigned as follows: 


ISTATE Processing: 
1 at the beginning of a Field. 
2 after the beginning but before the 


decimal point. 
3 after the decimal point. 
after the exponent sign. 
5 after an error has occurred in the 
field. 
Again, for different applications, the programmer may define 
other states as needed. 

The aext step in processing a character from the 
card-imag2 array is two-fold. The proper action to assemble 
the variable must be executed. This, of course, is given by 
the value of IACT. Also, if required, the state must be 
changed. & single integer can be used to represent both the 
next stats andthe assembly action required. Por example, 
the inforration contained in ISTATE and IACT can be combined 


into a single variable, IPROC, by: 
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EPROC = GACT* 100 te ES DATE 
The inforaation can be separated by: 
IACT = IPROC/100 
and 
ISTATE = MOD (IPROZ,100) 
A natrix can be assembled that has the row index defined by 
ISTATE and the column index definei by NSYM. If an element 
9f this matrix represants IPROC for the ISTATE and NSYM 
that define the element, then this matrix is an efficient 
method of programming the assembly task. The matrix so 
constructed is called tne Symbol-State Matrix. Reference 
(4) contaios futher discussion of Symbol-State matrices and 
their uses. 
For illustration, consider the problem of reading a 
real linear array of length NEL. Por Simplicity, elements 
of the array may not have exponents. The state iS 


identifiei as: 


LS PALE Processing: 

1 at the beginning of the field. 
2 before the decimal point. 

3 after the decimal point. 

4 after an error. 


The steps that ne2dad to be programmed are: 

a) [nitialize the storage array, the syaboli-state 
ec Gly, 2 hollerith translation table and the array 
pointers. A hollerith translation table is a vector that 
Contains the hollerith character indicated by NSYM in the 
element with the index NSYM. 

b) Read a card in Al format into a temporary integer 
array. 

c) ITiitialize the card column pointer called ICARD, the 
program state identifier (ISTATE), and the fractional digit 
counter called NCOUNT. 
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d) Iientify the character occupying the ICARD column of 
the card-image array (NSYM) from the hollerith translation 
table. 

e) Using NSYM and the current ISTATE as indices for 
entry in the Symbol-State matrix, called ISST, obtain IPROC 
Meer OC = [SST (ISTATE,NSYM). From IPROC, compute [TACT and 
the new ISTATE. 

f) Carry out the appropriate assembly action based on 
the value of IACT as follows: 

1) Set a flag indicating the number is negative. 

2) Assemble the digit into the integral part of 
the number (X = 10.0 * X + DIGIT). 

3) Assemble the digit ynto the fractional par tof 
the number (X = X + DIGIT/(10.0**NCOUNT)), and step the 
fractional digit count2r (NCOUNT). 

4) Print an error message. An action to 
compensate for the error may be included here. 

5) Store the element in the storage array. If the 
meray is filled, return; otherwise, reset NCOUNT and 
continue. If the symbol causing this action is EOR, reset 
ICARD and continue at step B. 

Note that a plus sign or d2cimal point changes’ the 
State of the program but does not require an assembly 
action. After the assembly action above is complete, proceed 
to step 6. 

6) Increment ICARD and test for the logical end of 
record. If ICARD is greater than the record length, reset 
ICARD, execute step £.5 and continue at step b. If ICARD is 
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less than or equal to the record length, continue at step d. 


re. 
SYABOL COLUAS DIGIT ¢ = ° OP EOR OrHeas 
STATE ROW 1 2 3 4 5 6 7 
BOP /BOR 1 202 602 102 602 501 501 aQu 
DIGIT BEPORE DZCISAL 2 202 404 404 603 $01 501 804 
DIGIT APTES DECIMAL 3 303 404 4O4 $040 501 $01 808 
ERROR 4% 604 604 604 604 601 601 6048 


SYSBOL ~ STATE TABLE 
TABLE 2 


Table 2 is a Symbol-State matrix that the programmer 
Might use to accomplish the variable assembly described 
above. The table presunes that the action to compensate for 
ah error is to set the element equal to 0.0 (or some other 
value) and execute step 5. Action 7 is assumed to be the 
same as 5 except it allows for differentiating between an 


end of record symbol and an end of file symbol with a test 





Sn the valaoe of IACT. This table does not indicate an error 
if there i5 no decimal point. Instead, the decimal point is 
assumed after the last digit of the field. Additions to the 
Symbol-State Table or modifications to the logic are not 
difficult once the concept is understood. 

Now consider the problem of assembling a character 
String, for instance, as a matrix name. If there are 4 
bytes per integer word, then one element of the card-inage 


array may contain: 


So) Ee ae a | 
where "bl" represents the blank character. The byte 


containing the hollerith charact2r S from the card-image 
array must be moved into the appropriate byte of the 


character string variable being assembled. This is nota 
trivial problen. ANSI FORTRAN does not support byte 
manipulation. Major computer aanufacturers have added 


features to their FORTRAN which allow byte manipulation for 
string assembly. CDC allows for byte manipulation with MASK 
and SHIFT operations. In IBM FORTRAN the job can be done 
with a suoroutine as follows: 


SUBROUTINE SHIPT(A,B,I,J) 
LOGICAL® 1 a (1) Bi1 
c----- THIS SUBROUTINE MOVES THE JTH BYTE OF 
= VARTABLE B INTO THE ITH BYTE OF 
: EARIAGLE A. 
RETURN 
END 


The next update of the ANSI FORTRAN Standard tentatively 
includes variable word size specification with a 
"CHARACTE2 *n" declaration, where nis an integer. If the 
change materializes then the abov2 subroutine will work for 
any pachine by changing LOSICAL*1 to CHARACTER*1. 
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C. DATA MANAGEMENT 


The comments on data managemeat have been restricted to 
the management of data storage in core. Only one of several 
alternatives is discussed and the management systen 
presented is not necessarily the most efficient method. It 
is advocated because it represents a Straight forward 
approach and one does not have to be a computer scientist to 
implement the procedure. 

The problem of in core data aanagement can be reduced 
to the problem of storing, finding and deleting arrays in 
the main array mentioned previously. One way of organizing 
the main array is to store an'array directory at the start 
of each arcay in the main array. The array directory needs 
Mmemcontaia the following information: 

2) The number of elements in the array. 

b) The precision of the array (integer 
words/felenent). 

s) A unigue identification for the array. 

d) The number of rows and/or the number of columns 


in the array. 


Diagram of Array Storage 


The process of reserving storage for an array should: 
3) Ensure adequate space iS available in the main 
arcay. Processing must be interrupted for resolution of the 
problem if sufficient room is not available. 
b) Adjust NSIZE to comp2nsate for the additional 
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wena y. 

cs) Create and store the array directory based on 
information from the calling program. 

d) Check boundary alignment and adjust NEL for the 
previous acray if necessary. It is possible to adjust NEL of 
tha curreit array so that any subsequent array will pe ona 
proper boundary no matter what the precision of the 
Subsequent array. This is a Simple solution but does make 
some available space unavailable. 

e) Return the main array address for the first 
eleament of the array tos the calling progran. 

Locating an array requires searching the array 
Mimeectori=s for awwatch vith the identifier of the array 
One methoi of doing this type of matching is: 

Te (cOmp (L(N) ,DE)) so TO nRnh 
where 
N is the address, in the main array, of the first 
word of the array identifier being tested. 
CI contains the name to be matched. 
ano is the statement number that is branched to if 


4 match 13s found. 


STOMP is a function subprogram defined as follows: 

FUNCTION COMP (I qh 

DIMENSION I (3) ,J (3) 

LOGICAL COMP 

COMP = .FALSE. 

DO 10 K = 1,3 

tine Sey ~NE. J(K)) RETURN 

. 10 CONTINU 

COMP = .TRUE. 

RETURN 

END 


This example assumes that three integer words are used to 
Sontain the array identifier. [If the arpay is found, the 
main array address of the starting element and the directory 
mufermatic: needs to be returned to the calling programs If 


the array is not found, processing should be interrupted for 








mesolution of the error. 

The motivation for providing a method of deleting arrays 
is two-foli. If matrices are not stored with unique names, 
th2 locating system described abov2 will only find the first 
array in storage whether or not it is the array needed. By 
uSing the delete op2ration under program (vice user) 
control, it is possible to ensure each array name is unique 
at the tim? it is created. The second reason for having a 
delete ciaoability is that when an array has served its 
purpose, deletion unier either program or user control 
allows the storage to be freed for other uSes. 

It 153 advantageous to leave all free storage postions 
at the eni of the alocated storage. The simplest way to do 
this 1s to identify, in the main array, the address of the 
first position used by the array to be deleted (N1), and the 
ficst elenent of the n2xt array in storage (N2), then: 


N3 = N2 - N1 

DO 10 I = N1,NSIZE 
10 L{I) = L(T#N3) 

NSIZE = NSIZE - N3 


Note that if the directory is stored preceeding the array, 
then N1 would refer to the first element of the directory. 
In tnis case an alternate way of computing N3 is: 
N3 = NEL*IPR + NDIR 

There is no reason to interrupt processing for error 
resolution when an array is not found to be deleted. It is 
helpful to print a message when an array is actually deleted 
especially if the DELEIE operation is initiated by the 
program, for instance, to prevent duplicate array names. 

The more complicated storage schemes for banded or 
skyline storage of arrays is also possible with this systen. 
Banded ani 

Skyline storage schemes make more efficient use of 


Storage by simply not storing many of the matrix elements 
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that are zero. The skyline storage algorithm is the most 
storage efficient of the two types. Reference (1) contains 
a substantive discussion of the skyline storage method, and 
examples of matrix manipulation subroutines for use with the 
Skyline algorithm. Suffice it here to say that the skyline 
storage alogrithm reguires the creation of two vectors. 
On2 is a vector of the elements from the matrix that are 
actually stored. The second vector contains the indices in 
the first vector of the diagonal elements of the matrix. 
Por this case, it is convenient to use the same directory 
for both vectors. If there are NEL elements under the 
skyline, 2ach element takes IPR integer words and there are 
NR rows in the matrix, then the storage required is 
NDIR +¢ NR + NEL * IPR 
a typical block of the main array might look like 


DIRECTORY | DIAGONAL ELEMENT INDICES STORED ELEMENTS 


The above presents a Straightforward approach to data 
Management in core. It is codeabl2= in ANSI PORTRAN IV by an 
lodividual familiar with the basics of the language. rt 
does not require a professional programmer for 


implementation. 


D. LOOPING 


The capability to repetitively execute a series of 
instructions (looping) is an important tool in problem 
solution. Frequently, for non-linear problems, there is no 
direct method of solution. An iterative scheme is normally 
employed to solve this type of problem. Brute force 


iteration is possible with any program by making multiple 
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ruas of the program with suitable alteration to the input 
data for 2ach run. This is very wasteful of both the user's 
time and computer time. On the other hand, if looping is 
possible, the user can program an iteration scheme using the 
commands available, then place the series of commands in a 
loop. A relatively simple method of providing a looping 
Capability is to store the commanis of the loop as a matrix 
in the main array. Then, with appropriate pointers, it is 
possible to execute a series of inStructions repetitively 
from memory, just as if they were coming from the input 
device. Details of the method of looping used in the 
program CAL are contained in Chaptar IIT. 

A particularly attractive advantage of looping is the 
ability to iterate for selected eigenvalues of a matrix. 
This ability permits one to ‘approximate the number of 
Significant digits in the computer soluton of a problem as 
follows: | 

NSD = NDA - ALOGIO(EMAX/EMIN) 
whera 

NSD is the number of significant digits. 

NDA is the number of digits available in the word 
size. 

EMAX is the maximum eigenvalue of the matrix. 

EMIN is the minimum eigenvalue. 

ALOG10 is the ANSI standard function for the 
logarithm to the base 10. 
This approximation holds for positive, definite, Symmetric 
matrices. Computer results for problems where NSD 


approacnes 0 should be viewed with extreme skepticism. 





Tit. CAL ORGANIZADION 


This chapter provides a general overview of the 
operation and organization of CAL. It is intended to 
provide sufficient information td permit users to progran 
Subroutines for use in CAL and to operate the CAL systen 
sensibly. 

CAL iS a modern computer program written in the FORTRAN 
IV programming language. It was written as a teaching aid 
for illustrating the direct stiffness method of structural 
analysis. in addrt:on to extensive linear algebra 
capabilities, the program provides several analytical 
alternatives for both static and dynamic analysis of elastic 
Structures. 

CAL execution is flexible, controlled by user selection 
of operations in a logical sequence (not unique) from the 
Operations available in the prograna. fhe input data deck is 
a seguence of modules. The first card of each module 
contains the name of the requested operation, a list of 
Matrix aryument names, and integ2r parameters used in that 
Operation. Comments may be placed on the card following a 
blank. 

The »9peration cardi looks like this: 

TPP Rar rON At, Zee. NI, Ninn NE | =6COMMENTS 
M1 throug: M9 are matrix argument names. A matrix argument 
Name consists of up to 8 alpha-numeric characters, the first 
being alphabetical. N1 throuch N4 are positive integers. 
An operation card may use from zero through nine matrix 
arguments and/or from zero through four integers. If 
required, the operation card is followed by a sequence of 
data cards containing numerical information used for the 


step. Some features of the output from CAL are aiso 
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controlled by user selectable operations. Details of the 
Dperations available and how to execute CAL at NPS are 
presented in Appendix B. 

CAL 15 logically iivided into four segments. The main 
segment contains the main program and service subroutines 
used comndonally by all other s2gments. In subseguent 
descriptions, the main segment will be referred to as CAL. 
Among other things CAL contains the subroutines to deliver 
the next operation and the subroutines for data management. 
The three cemaining segments will be referred to as GROUP 1, 
GRIUP 2, and GROUP 3. 

GROUP 1 contains subroutines which perform matrix 
mimepalations, i.@, load, print, -aWltiply,. etc. GROUP 2 
contains subroutines associated with formulation of the 
Static problem and display of the computation results. For 
example, there are subroutines to: . 

a) input nodal geometry; 

b) input boundary conditions, 

=e input Loading conditions, 

41) form element stifness and amass matrices by 
several different methods; 

2) combine element matrices into global matrices; 

f) display results. 
Finally, GROUP 3 contains subroutines associated with 
dynamic analysis. Two methods of evaluating the second 
order egiations of motion, eigen analysis, and printer 
plotting routines are available. 

Matrix storage in CAL is don2 dynamically as described 
in Chapter II. It is not possible to state with precision 
the maximum problem size. Experience shows that with the 
default parameters of the NPS version, static analysis 
problems ace limited to approximately 48 degrees of freedom 


(DJF). Dynamic analysis problem size is controlled by both 
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the system DOF and the number of jlisplacement vectors to be 
calculatei. Por a 25 DOF system, about 100 displacement 
vectors caa be calculated. 

All acrayS and matrices in CAL are stored ina one 
dimensional matrix specified in blank common as L(XXXXX). 
XXXXX represents the total storag2 reserved for subscripted 
Variables. Each array is preceeded by a directory which 
contains: 

1) the number of elements in the array (NEL). 

2) the number of integer words required to store 
an element (IPR), 

3) the number of columns in the array, 

4) the number of rows in the array, 

>) the first four characters of the matrix name, 

and 5) the last four characters of the matrix name. 

An array uses (NEL*IPR + 6) integer words of storage. 
Larger problems can be solved by increasing the size of the 
L array in the main program. Th2 variable MAX must be set 
to the surrent dimension of the L array. This data 
Management system is described in Chapter II. CAL uses 
three subroutines for the data management sub-system. There 
1s a subroutine called LIST which is used to reserve storage 
for new matrices as they are creat2d. The calling arguments 
are, in orier, NM, NR, NC, IPR, and IERR. NM iS’ an array 
mane assigned by the user. Array names must start with an 
alohabetic character and conatin 1 to 8 alpha-numeric 
Characters. 

Since the name, NM, is us2i only in the directory, 
there is no reason to conform to FORTRAN name conventions. 
Thus, names beginning with the characters I through N may be 
asSigned for real matrices with impunity. Call the 
Stiffness matrix K, the mass matrix MASS andthe nodal 


coordinat2 array NODES. All variables are treated as real, 


36 


Oe ee ee ee ee ee | = = 98 8 ee ow ow 








ee oe wage ow ~~ 


double precision humbers unless clearly specified as 
integers in the command description. 

NR and NC are respectively, the number of rows and the 
number of solumns for the matrix to be stored. IPR is the 
number of integer words required to store an element. IERR 
1S an error parameter. iferekheas not daqual to 1, then»an 
error has occurred during the list operation. 

The subroutine LOZTATE is usei to locate arrays in the L 
arcay. Th2 calling arguments are, in order, NM, NA, WR, NC, 
ani IERR. WM, NR, NC, and IERR are used as described above. 
NA is the index in the L array for the first element of the 
arcay. I1 a subroutine call, L(NA) is used to pass the 
matrix NM to and from the called subroutine. 

The remaining data management subroutine is DELETE. It 
has asingle calling argument. NM 1S used as discussed 
previously. DELETE removes the matrix NM from the L array 
and frees the storage taken by NM at the end of the L array. 
Figure 7 illustrates this process. 


Directories D, and Matrices M 






DELETE M 
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One Dimensional Array Storage 
FIGURE 7 


NOw consider once more the function of the main progran 
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segment CiL. In addition to the data management subroutines 
described above, the CAL segment contains four subroutines 
us2d to deliver the next operation to the calling progran 
segment, anda timing subroutine. The subroutine INPUT 
determines whether the next operation is part of a loop or a 
new command from the input device. If the next operation is 
not part sf a loop, INPUI calls the subroutine OPREAD. In 
tucn, OPREAD calls the subroutine RCARD to read and 
interpret the free format inforsation of the next command. 
Pree format data input is discussed in Chapter II. [If the 
command initiates looping, then OPREAD calls RCARD until all 
operations in the loop have been read. As operations are 
received by OPREAD, this subroutine stores them row-wise in 
the L array. After the last operation of the loop has been 
read, OPRZAD calls LIST to prepare the directory. INPUT 
calls JPSfIR to provide each command in seguence from the 
loop matrix inthe L array until the loop operations are 
complete. The subroutine SECOND is called by INPUT to 
determine the amount of time used in each computation step. 
A block diagram of the connections is: 





SECOND 





OPREAD 








Cronk 


A bl»ck diagram of the total system follows: 
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MAIN INPUT 


sROUP1 GROUP2 GROUP 3 
INPIJT INPUT INPUT 
ANALYSIS ANALYSIS ANALYSIS 


From this example, it is apparant that the program is 
Organized so that only one group at a time needs to be 
resident in the core with CAL during execution. Thus, an 
Overlay structure can be used which permits CAL to operate 
in a very small computer system. 

Although CAL can operate ina small (64K) computer 
System, a fzass storage device (magnetic tape or disk) must 
be available both to store program overlays and for use as a 
scratch pad during structural analysis operations. Element 
stiffness and mass matrices are written out on a amass 
storage device as they are created, then read back into the 
computer nemory as needed for assembly into the global 
matrices. This method minimizes the amount of storage 
reyuired in the core. Space is required to store the 
matrices only an element at a time, no matter how many 
elements are used. 

The NPS version of CAL also contains commands which 
permit a user to stop execution before the end of a problem 
then resume at a later time. The operations also use a mass 
storage device to read or write the contents of the L array 
and certain other parameters. 

One of the features of CAL which makes it so flexible 


is the availability of looping. Loops can be nested up to 


Sy, 














five deed. Loops are initiated by the LOOP command and 
terminated by the NEXT command. Por operations within a 
loop whith reguire data, the data cards must appear in 
proper sejyuence, but after the last NEXT card of the loop. 
An example is shown with the loop command description in 
Appendix 8B. A skip command is available which allows 
selective skipping of operations within a loop. 

A second feature of CAL whith qives it really great 
flexibility is that user coded subroutines can be called 
under projycam control. GROUP 1 contains two operation names 
reserved for user coded subroutines. These names are USERA 
and USERB. A user soded subroutine might beused, for 
instance, to compute equivalent nodal forces for a 
MeStrioutei load. 

Two subroutines are needed’ for a user operation. The 
actual numerical operation is coded. normally in a subroutine 
with any name not already used in the CAL program. A buffer 
Subroutin2? must then be written to interface the numerical 
computation with the CAL data management system. lle 
illustrate the technigue of interface, consider the familiar 
example of matrix multiplication (already available in CAL). 
Symbolically the problem is to compute C where C = A*B. The 
basic algorithm for matrix muitiplication is 

C(L,J) =  B&(T.K)*B (K, J) 
A subroutine for the job might be: 

SUBROUTINE MULTIPLY (A,B,C, NRA,NCA,NCB) 
DIMENSION A(NRA,NCA) ,B( NCA, NCB) ,C (NRA, NCB) 
DD 20 I=1,NRA 
DO 20 J=1,NCA 
KX = 0.0 
DI 10 K=1,NCA 

10 X = X + A(I,K) *B(K,J) 

20 ={I,J) = x 








RETURN 

END 
The opera:ion card might be selected as: 

USERA,M1, 82,43 COMMENTS 

where M1, M2, and M3 are respectively the matrix argument 
manes for the dummy arguments A, 8B, and C. The recommended 
st2ps to be programmed in the buffer subroutine named USERA 
Ee 

a) locate 41 and M2; 

b) ensure thea number of columns in A is egual to 
the number of rows in B; 

cs) call the DELETE subroutine to ensure that a 
matrix with the same name as M3 does not already exist in 
storage; 

d) call LIST to reserve storage for M3; 

2) test IERR to ensure that no errors have 
occurred in the above steps; 

and fy call subroutine MULTIPLY to complete the 

operation; 
At the minimum, steps a, d, and f must be programmed. Here 
1s a possible subroutine to meet the requirement. 

SUBROUTINE USERA (IERR) 

COMMON MAX,NDP,L(1) 

COMMON /CARD/ INHOL (3,10) ,N (4) 

CALL LOCATE (INHOL (1,2) ,IA,NRA,NCA,IERR) 

CALL LOCATE(INHOL (1,3) ,1B, NRB,NCB,1IERR) 

IF (NRA.NE.NCA) IERR= 2 

CALL DELETE (INHOL (1,4) ) 

CALL LIST(INHOL (1,4) ,NRA,NCB,NDP,IERR) 

CALL LOCATE(INHOL (1,4),1IC,NRA,NCB,IERR) 

IP (TERR.NE.1) RETURN 

CALL MULTIPLY (L(IA),L(IB) ,L(IC) ,NRA,NCA,NCB) 

RETURN 


41 





END 

LOCATE, LEST, and DELETE are used as described previously. 
Blank common appears sd that USERA has access to the L array 
and the parameter ND?. The variable NDP is initialized at 
2, the nuaber of integer words r2quired to store a double 
precision variable. The named common CARD is initialized by 
RCARD or OPSTOR. Th2 INHOL array contains the operation 
name and the matrix argument nam2s. Fach name is stored 
columnwis2, i.e., column 1 has th2 operation name (USERA in 
this case), column 2 has the name of thematrix argument M1, 
etc. Row 1 elements have the first four characters, and row 
2 elements contain the remaining four characters of the 
respective argument. Row 3 is not used in the NPS version 
SeCALbwe sho the farst call to subcoutine LOCATE, INHOL(1,2) 
provides the name assigned by th2 user to matrix argument 
ay. On@return €om LICATE, I* will be the index in the L 
arcay of the first element of M1. IB and IC respectively 
provide th2 same inforaation for M2 and M3. Other variables 
are used as discussed previously. 

Now, generalizing the proceiure of the above example, 
th2 user must program two subroutines to define a user 
Operation for CAL. One subroutin2 performs the operation, 
and the second acts asa buffer between CAL and the 
operation. Iwo names are available, USERA and USERB, for 
the buffer subroutine. Both have a single passing argument, 
the parameter IERR. The program assumes an error has 
occurred when IERR takes on a value other than 1. The 
buffer suyroutine must: 

a) locate matrix arguments designated for the 
Operation, 

b) reserve storage for arrays created by the 
Operation, 


Sc) pass the elements of the L array for the matrix 


42 





parameters used in the operation. 

In addition, it is recommended that the buffer subroutine: 
a4) ensure that the matrix arguments to be used are 

compatible with the operation, 1.¢., square, symmetrical, 

etc. 

b) call delete with names of matrices to  »b2 
created prior to reserving storage. This will prevent 
Matrices with duplicate names. 

c) test TERR to ensure n> error has occurred prior 
to calling the operation subroutine. 

Mas Operacion subroutine can be coded with normal PORTRAN 
techniques. 

In this chapter, the internal organization of CAL has 
been described, along with several important features of the 
Baogranz. Most o£ the non-numerical procedures described in 
the preceeding chapter are illustrated in CAL. This chapter 
has highlighted the flexibility of CAL and provided an 
example td assist users in coding their own operations to be 
executed by CAL. 
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ieee bene Nt ONG CAL ON RODIFFEREET COMPUTER 


In this chapter, the procedure to convert CAL for use 
on a different computer system is discussed. The NPS 
version of CAL has been organized to minimize the difficulty 
associated with implementing the program on a new computer. 
Th2re renains, however, some things that are machine 
dependent. 

The following information is required for the new 
Systen. 

a3) How many characters can be stored in one 
integer word? 

b) How many integer words does it take to store a 
real variable? For structural analysis, a real variable 
Should contain at least 12 significant figures. 

*s) What method can be us24 to manipulate bytes? 

1) How is elapsed CPU time computed? 

[Th2 answers to the above questions determine the complexity 


of the conversion task. 
We IBM FORIRAN STATEMENTS 


Ther are a few statements in this version of CAL 
solely in the IBM FORTRAN Language. The statement "IMPLICIT 
REAL*8 (A-H,O-Z)" makes all real variables double precesion. 
This statement appears in almost every subroutine. It 
Should be replaced by an appropriate statement when using 
another conputer. The calls to ERRSET suppress printout of 
error messages. These statements appear together in 
Subroutine CAL1 and must be removei. They do not need to be 
replaced. The call to SETIME in the main program should be 


renoved for use on non-IBM machines. 
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pee LIME COMPUTATION 


There 1S no standard method of computing elapsed CPU 
time. The subroutine SECOND must be rewritten to conpute 
elapsed C279 time in seconds for th2 new system. The library 
of subroutines for the computer system should have a 


Suitable routine. 


oe BYTE HANI PULATION 

As previously mentioned,’ there is no standard way to 
shift bytes in a variable. The subroutine MOVEB will work 
if one byte logical variables are available. If this is not 
th2 case, than MOVEB must be ‘rewritten using a byte 
manipulation technigue suitable for the computer systen. 


D. OPERALTION AND MATRIX NAMES 


In this version, names are stored four characters toa 
Variable since the IBM 360/67 computer uses four bytes per 
Standard integer word. Names are stored column-wise in 
vectors of three elements. The parameter NH must be set to 
the number of elements required for storage of a name on the 
new systen. LBYTE must be set td the number of characters 
per integer word. NH should be set such that NH*LBYTE at 
least 6. Both of the parameters are set in the BLOCK DATA 
Subprogram. Data statements initialize operation names in 
the BLOCK DATA subroutine andin all three of the GROUP 
subroutines. All these data statements must be adjusted so 
that they initialize LBYTE characters per element. Format 


Statements for write commands with output in format A4 
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Should be changed. For the most part the output format wAd 
should be changed so that w = NH and d= LBYTE. Other uses 
Should be apparent to the programmer. 


E. _BOUNDARY ALIGNMENT 


Ther2 are three potential areas which may cause 
boundary aiignment problems. NDP is set in the main progran 
to the number of integer words needed to store a real 
Variable. NDIR is the number of integer words used for a 
Maite Gax directory. The first possible problem arza is the 
relation between NDIR and NDP. If NDIR is not evenly 
divisable by NDP, add the remainder to the eguation for NDIR 
in the BLOZK DATA Subroutine. The second possible problena 
area is in the Subroutine LIST. his subroutine contains a 
test to ensure the number of elenents reserved for a single 
precision array is evenly divisibl2 by NDP. If NDP is more 
than 2, taz= test 1S not valid. A better test would be 

IP (N2.~EQ.%1) IT(1) = NR * NC + MOD(NR*¥NC, NDP). 
Pinally, if NDP is greater than 2, then blank common must be 
adjusted s> that there are NDP integer variables preceeding 
the L array in the main program and all subroutines where 


blank comzon appears. 
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Ve. CONCTLUSIONS AND RECOMMENDATIONS 


CAL provides a flexible tool for teaching modern 
sturctural analysis. The general matrix operations make it 
us2ful in the area of linear algebra as well as structural 


analysis. The use of CAL is highly encouraged. 


A. RECOMT ENDED MODIPICATIONS 


A graphical display capability should be incorporated in 
CAL as soon as the equipment .becomes available. It would be 
2spically helpful to have a display of the element forces 
from the FORCE operation. Additionally, incorporation of 
heat transfer options into the program should be considered. 
The heat transfer capability would be a valuable tool for 


Students in the department. 
Coro Lo ER BPPLICATIONS | 


The organization of CAL provides an excellent outline 
for development of other major analytical programs. The 
Storage scheme can be simply mddified to permit matrix 
storage using the skyline algorithno. Large problems can 
then be ascomodated by incorporating existing algorithms for 


an out-of-core equation solver. 








APPENDIX A - PROGRAM LISTING 
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This appendix provides details on use of CAL with the 
IBM 360/67 computer at NPS. Sestion I provides details on 
the command structure. Section II is a Summary of commands 
available. Section III provides the job control language 
for executing the program in both the batch and interactive 
modes at NPS. Section IV contains detailed specifications 
for each available command. Finally, section V gives 
direction for solving larger problems with CAL. The 
majority of this appendix. was originally published as 
rezterence (3). The author wishes to express appreciation to 


Professor Wilson for permission to use this material. 
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foe LORMAN D RESTRICTION OF THE LANSUAGE 

CAL is an interpretive language which is designed to 
Manipulate arrays and matrices and to perform several 
Standard structural analysis operations. A CAL program run 
involves the reading of the input deck once and executing 
the commands designated by the operation cards as they are 
2ncountered. Loopiag Operations allow a seguence of 
commands tS be executel more than once. 

The input deck is composed of operation cards and data 
cacds. The data cards directly follow each operation card 
which reguires data (see LOOP operation for exception to 
md S) . [he operation card -:contains the name of the 
Operation t2 be executed, names of arrays associated with 
the operation and integer constants. Examples of the 
general form of this card are | 

OP, M1,M2,43,44,M5,N1,N2,N3, NG COMMENTS 

OP, M1,N1,N2 

OP,N1 

oF 
in which JP is the name of the op2ration to be executed, Mi 
is the nam2 of an array and Niis an integer. The names of 
9P or Mi are one to eight alphabetic or numeric characters 
to be selected by the user. The first character of a name 
must be alphabetic. The sequence of terms OP,Mi and Ni must 
be separated by commas. Characters following a blank will 
be printeil aS comments in the output from the program run. 

If an Operation attempts to load or generate an array 

which previously existed the program will delete the array 
before the execution of the operation. A new array need not 


be the same size of the old array which had the same name. 
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II. 


A. GENERAL COMMANDS 


fernadicdt2=s a Significant change ofr addition in CAL=-NPS 


START 
STOP 
NO 
LES 
LABEL 
READ 
WRITE 
TIME 
SAVr 


2 Je by 


ae > T 


Initialize for the next problen 
Normal termination 

Temporary supprossion of sutput 
zestores output 

Print comments | 

shange logical device for input 
Shange logical device for output 
Suppress time printout — 
Interrupt a problen . 
>Ontinue an interrupted problen 
List arrays and storage used 


B. GENERAL MATRIX COMMANDS 


Load user defined matrix 
create null or unit matrix 
Yatrix print operation 

Matrix duplication 

Yatrix addition 

Matrix subtraction 

Matrix multiplication 

Yatrix transpose 

“ultiply a matrix by a scalar 
Solution of linear equations 
Form sub-matrix from large matrix 
Store sub-matrix in large matrix 


Form row matrix from diagonal 


Jane 








STIDS 
MAX 
NORM 
INVEL 
SQREL 
LIS 
PROD 


DELETE 


store row on diagonal 

Evaluate row maximums 

Evaluate matrix norms 

Inverte each term in matrix 

Sjuare root of each term in matrix 
Natural log of each term in matrix 
Evaluate product of all terms in matrix 


Delete matrix from storage 


ce oTATITC ANALYSIS QJPERATIONS 


NODES 
BOUND 
BEAM 
2RuU SS 
PLANE 
SLIPE 
FRA ME 
LOAD 
LOADS 
ADDSF 
ADDK 


MEMFRC 


DIS PL 
PIRCE 


D. 


FUNG 
Sit ey 
EIGEN 
DYNAM 


PLOT 


Input joint geometry 

Specify boundary conditions 
Porm 3-D beam stiffness matrix 
Porm 3-D truss stiffness matrix 
Porm 3 to 8 node plane stiffness matrix 
79rom stiffness matrix froa slope/deflection eq. 
?orm 2-D frame stiffness matrix 


Load integer array (connectivity matrix) 


Form load vector 
Porm global stiffness and mass matrices 

idd element matrix to global matrix 

Calculate element forces from joint displacements 
Point joint displacements 


Evaluate and print member forces 


DYNAMNC SY ANALYSTS OPERATIONS 


- 32nerate equal interval time function 


- Integrate dynamic equilibrium equations 


2valuate mode shapes and frequencies 


- Evaluate uncoupled e eons of motion by 


-= 


mode superposition metho 


Line printer plot of joint time history 
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E. LOOP JPERATIONS 


EOOP =e aGt Of LOOD 

NEXT Zoid NOLt Loop 

SKIP svsogdaitional skip of operations within loop 
E NAMES AVAILABLE FOR USER SUBROUTINES 

OSERA 

USERB 
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ft. JCL FOR EXECUTION 
The jbb control language (J<OL) necessary to initiate 
program execution for the OS/MVT (Datch) system is: 


//({standard green job card) 
f/f EXEC CAL 
J/SYSIN DD * 


CAL COMMANDS AND DATA CARDS 


/* 


For the CP/IOMS (time-sharing) system the procedure is: 


use standard login procedure, then 
cP LINK OQO99P 191 199 PASS= ESAN 
LOGIN 199 A,P 

CAL 


CAL COMMANDS AND DATA CARDS 


The following information is provided for the user who 
wishes to use his own subroutine or solve larger problems 
then allowed with the default core allocation. The load 
module library for the OS/MVYI operating system is in 
partitioned data set F0099.CAL on DISKO2. The load module 
is assigned member name MODULE. [he file used as a scratch 
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Bad 1S FrJ1P001. Pile FPTO2F001 is used as the scratch file 
for the SAVE and RESUME operations. The PTO1F001 file is 
currently pointed at the SYSDA disk and FTO2F001 is pointed 
at dummy. If for any reason a us2r desires to use the SAVE 
Operation with the IS/MVI system, both FTOQ1P001 and FTO2FO001 
must be redifined. Ses Computer Canter TN No. 0141-05, USER 
LIBRARIES AND SOURCE CODE EDITING UNDER OS for procedures. 
File FIO1F001 requices 3404 bytes of storage per element. 
File PTO2P001 requires a total of 4 times the dimension of 
the L array plus 8 bytes and should have RECFM=VBS. 
For C2/2OMS, the execution routine is 
ETYPEQUT OFF : 
CP Sst LINELN 130 
GLOBAL T SYSLIB SSPLIB 
PILEJEF 01 DSK-P4 FILE FTO1F001 RECFM VS LRECL 3408 
BLKSIZE 3412 2 
FILEDEF O02 ODSK-PU FILE FTO2F001 RECFM VS LRECL 7290 
BLKSIZE 7298 
Poeeceoge Ou DSK 
FILEDEF O08 PRT RECFM FPBA LRECL 133 (PERM) 
LOAD STORE USER CALH (CLEAR NOMAP XEQ) 


FTQ1P001 15 the scratch pad for the structural analysis. 
FTO2FP001 15 created by the SAVE operation. Both FTO1TFOO1 
and FTO2F)01 must be available for the RESUME operation. 
Temporary additional storage can be obtained, if needed, 

by typinjy GETTEMP while in CMS. The GETTEMP execution 
routine is: 

&TYPEOUT OFF 

VSET ROYSSG OFF 

CP Sil LINELN 130 

&SBEGSTACK 

Worl BLIP + 
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CP DERPINE T2314 192-4 
FORMAT T ALL (NOTYPE) 
RELEASE 192 T 

LOGIN 192 P 

LOGIN 191 B,°? 
GENDSTACK 


FOS 








IV. CAL COMMAND SPECIFICATIONS 
A. GENERAL MATRIX OPERATIONS 


CAL has most of the standard matrix operations plus 
some special array operations which are useful in structural 
analysis. The following is a list of approximately 25 
Op2rations which are used for control and general matrix 
manipulation. 

+ indicates the formation of a new matrix. A matrix 
previously defined with the same name will be deleted. 


- indicat2s modification of an existing matrix. 


_—eap ween @weaweamaeDbee ee ee SPs es SFP SF SP FSF PF SPP ZP as BSPeese ses SF SF 2 OS FP BZ Vs SF SPs Se BF SS SP Ss Zs Ss S&S = = 


This operation eliminates all arrays which were perviously 
loaded or generated. 


reese FFP eee eax2e es Fae FSP at BP SPP Se SP SPO Se Se SP SP PSP SS ee SS ee eee ee See ee eee ee ee ee See oe ee © as ae 


Th2se op2rations are used to selectively suppress output 
from CAL. The NO operations suppresses all printing, except 
Gesqnosticscs, wnatil the operation YES 1s encountered. 
Therefore, in subseguent runs of the same CAL progran 

Sutput which was previously correct need not be reprinted if 
these cards are inserted in the data deck. 


LABEL,N1 


This operation will read and print Ni comment cards, which 
follow this operation card. Colunn 1 of each card will be 
MimeeapEecteCdnas a Standard™carriag= control sygbol (i.e. 0 
for doubl= space and 1 for skip to the top of the next 


page). 
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= eS aaa IS VALID ONLY WITH THE CP/CMS TIME SHARING 
This operation permits the selection of the offline printer 
or the terminal as the input file device. Default is the 
terminal. N1 =4 will read subseguent commands from FILE 
FTO4FOO1T on the. cede Sr ee N = will restore the 
terminal as the input file device. All disk files prepared 
for use w#ith this command should end with either STOP or 
READ,4. [his command will not be executed onthe OS/SVT 


qwe@wmpeeeee SB wee BP ee ese es BPS ee ass eet ese 8S eS Se BPP eee EES Sse ee ee See ee Se ee ee ee ee ee ee ee 2 2 a 


WRITE, Nt 
ecco WS VALED ONLY WITH THE CP/CA&S TLIAE SHARING 


This operation permits the selection of the offline printer 
or the terminal as the outpyt file device. Default is the 
terminal and all error messages will be printed at the 
terminal Bede tess of the output file device selected. Ni = 
3 selects the offline printer. N11 = 6 restores the terminal 
as the output file device. This command will not be 
executed on the OS/MVYT (batch) systen. 


TIME | 
This operation permits the time printout to be suppressed 
without 1953s of other output. |. A second TIME will restore 


the time printout unless the print output is suppressed with 
the NO command. 


The SAVE command creates PILE FPTO2F001 containing all arrays 
lo storag2, at the time, of issuance. Note that Saved arrays 
Will contain any modification since creation. For instance, 
1f a matrix has been reduced by the SOLVE operation, the 
reduced forn of the matrix will be stowed in FILE PTO2P001. 
Note the reguirements discussed previously for use of this 
command with the OS systen. 


mepanp eee eee eee we ee @ 2 Pape a em one ee 2 SP ae 2 ae eae a 2 ee 2 eee | OP ew ew SE eee ee SS a aw SE SD SP Se Ses 2 Se = C= aD 


RESUME 
The RESUME operation reads FILE FIO2F001 into memory. An 
arrays cucrently in storage will be destroyed. Fit 


FTO2P001 must have previously been created on a _ mass storage 
device using the SAVE operation. Note that FILE FTO1TFOO1 
must also be accessible if an interrupted structural 
analysis problem is being resumed. 


Th> LIST dXperation prints the jirectory information for 
arrcayS in storage and the amount of storage used. 
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+ 
LOAD, M1,N1,N2,N3 


This operation will load an array of real numbers named 41 
which has Ni rows and N2 columns. The terms of the array 
ar2 punch2i in row-wise sequence on data cards following 


this operation. If N3 is_ zero or blank, the cards are 
punched ina format of (8F10.0). Tf @3 1S nonzero,” an 
additional card containing the format of the data cards must 
follow this operation and preceed the data, cards. If the 
data 1s to be numbers per card in field widths of 15, the 
eis tone cards would contain the following information: 


+ 
ZERO,M1,N1,N2,N3, NG 


A real natrix named M1 is created. with N1 rows and N2 
columns. [he terms in this-matrix will have the following 


values. 

Bey (1, 1) N3 SAS So. 

Bla, NG J = Iye6- : 
Th2refore,: this operation can be used to form null or unit 
matrices. 


=~ aes eee Fetter SS SP SS SBP eee zagse ewe ees B22 SP PP FS BF ese BPP SPS SPP SP eewaew SS SS SSS a Ss a oD 


Beonr,M1 at PRINT,M1,N1 or PRINT,N1,N1,N2 or PRINT, M1,,N2 


This. Se Nec will print. _ the r2al_ array named, M1in a 
Matrix format of up to 8 columns per line. If N1 is greated 
than zero the operation will read and print N1 comment cards 
wich follow the operation card. .N2. 1S optional but note 
that an extra comma nust be used. in place of N1if_ no 
printed title is desired. The matrix M1 will be printed in 

actitioned form with N2 columns per partition. Lines will 
ave N1*15+5 characters. N2 defaults to 8, peg ing 125 
characters per line. The user is cautioned not to overcome 
the capacity of the printing devic2 in use. 


ee “7 moor neo BA SS ee ee 
DUP, M1, M2 


This operation will foro 


an array named M2 which is 
l1dentical to the array naned M1 


—_ ee em em oe eo ee ee ee Ss ee ee ee eo a sa a es a ae ae ee ee i i is a si ais i iia oa ee ee ee 


ADD,M1, M2 


This operation will replace matrix M1 with the sum of the 
Matrices 41 and M2. 


SU0B,M1, M2 


This operation will replace matrix M1 with matrix M1 less 
Matera x 2. 
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+ 
MULT, M1,u42, 83 


This operation generates a new 


in) * M3 which is the 
product of matrices M1 and M2, or 4&3 1*42 


+ 


This operation generates a new matrix M2 which is the 
transpose of matrix 41. 


penoboe, it, 42 
This operation replaces each t2rn in the matrix named M1 


with the term multiplied by tha tarm M2(1,1) of the matrix 
naned M2. ‘ 


mE_ ee eeelUL6EweEUUlUlUeel eee eer UU eer Ue Ue Umer lle ee ee ee 626 ee am a 6a ae =e apnaee eee epee aS SB Se eS SS S22 SP SP HBSS SS Se SE Se ee ee Se ee ee an 4» == 


Seve tiie, Nien 2 OF SOLVE,M1,82,,N2 or SOLVE,M1,N1,N2 or 
SOLVE,M1,M2,N1 Om SOLVE, 41,42 or SOLVE, M1 

If Ni = 9), this operation solves the matrix equation AX=B. 
Mi is the name of the A matrix.,and M2 is the name of the B 
Matrix. Matrix A is trangularized. and the results X, are 
Stored in M2. . 

If N= Matrix A is triangularized only. 


1 
 Ni=2 For a given B matrix and the Matrix previously 
tciangularized, fhe B matrix is replaced by the results X. 
=3 Matrix A is replaced by its inverse FOR SYMMETRIX 
MATRICES ONLY. ; 
N2=0 or blank, matrix A.1s symmetric. If N2 1S nonzero 
the matrix A is not symmetric. 
For symmetric matrices, matrix A is factored intd the LDL 
form. The iiagonal D matrix is stored on the diagonal of A. 
The parameter N2 ermits the direct solution of 
non-symmetcric systems o equations. If N2 is not equal to 
0, an LU jiecomposition of matrix A will be performed. No 
direct replacement of 41 by its inverse is availabe for the 
non-symmetric case. Instead, use the ZERO operation to 
create an ea matrix M2 of the same order as Mi. The 
command SOLVE,M1,H2,,N2 will then replace the matrix N2 with 
the inverse of the matrix A. 


anpeanp 4) 42 == «== «=» —_— oe ee ee ee ee es eo ee a em oe Se a ee Ss ee Se SP eS ee ee 2 ee es ee SS ee ee a a a a a 


> 
DUPSM,M1,M2,N1,N2,N3,N4 


This operation forms a4 new submatrix named M2 with N3 rows 
and N4 columns from terms within the matrix named M1. The 
first tere of matrix A2, Be AS 1),. will be from row N11 and 
column N2 of matrix M1, or 4M (N1,N2) 
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sesh, 1,42,N1,N2 


This operation stores a submatrix named M2 within the matrix 
named 41. The first term of the Submatrix M2 will be stored 
at row N1 and column N2 of matrix M1. The terms within the 
area of Ml in which M2 is stored will be destroyed. 


& 
DUPDG,M1, M2 


This operation forms a new row matrix named M2 from the 
diagonal t2rms of matrix M1. 


ow ewe @& = ee Sauwenewreaeart se Se Be ease SPP SP Se SP SSeS SPSS ee Se ee SB SP SP ee ee eee ee ee 2 2 ee See Se eee Se Se ee ee a a =a» 


STIDG,M1,M%2 


Tnis operation stores a row or column matrix named M2 at the 
diagonal Locations of matrix Ml. 


+ 
MAX ,M1, M2 


This operition forms a column matrix. named M2 in which each 
row contains the maximum absolute value of the corresponding 
Mowe ch mitcix M1. The maximum and its column number is 
printed ror each row. 


+ 
NORM,M1,M2,N1 


If N} =), a, Tow matrix named M2 is formed in which each 
column contains the sum of the absolute values of the 
corres Bhp oule column Of Miatrix M1. Lf Ni /= 0, a@ crow gpatrix 
named M2 13 formed in which each column contains the square 
root of the sum of the squares of the values of the 
cocresponiing columns of matrix M1. 


INVEL, M1 


This operation replaces each term in the matrix named M1 
With it's inverse. 


SOREL, M1 


This operation replaces each tern in the matrix named M1 
with the sgyuare root of the term. 


nee eee wens ss ese 2 Se SF Se Se Rae Fe ee ese ee eSB SS SP ee ee ee Se eS Se SP ee ee eee ee ee See See ay Se = = oP == 


LOS ,M1 


This operition replaces each term in the matrix named 41 
With the noatural log of the tern. 
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+ 
PROD, M1,M2 


This operation forms a4 1.x 2 array named M2 which contains 
the product of all terms in the matrix M1. The product, Xx, 
is stored as two numbers of the forn 

J X = P*10**E 

moewouch 42(1) = P and M2(2) = E, the emponent. 


mPa e 8 ae SBS PF SPP eee SS eee ee ese BP ee ee ae BP aS SP eee ae Se See Se Se eee ee eee ee Se ee Se ee SP ee Se ie ee ee Se Se ee ae Se ee & ae 


PECETE, M1 


This operation will cause the elimination from storage of 
the array named 41. 
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Be. STATIC ANALYSIS OPERATIONS 


The purpose of this series of operations is to form-the 
total stiffness and diagonal (lumped) mass matrices for 
systems Sf two eng three-dimensional elements. For 
three-edimensional analysis, there are beam and truss 
elements available. For two-dimensional analysis, there is 
a frame elament, a slope/deflection element for beams, and a 
3 to 8 node isoparametric finite element available. 

After the creation of an array containing the 
coordinates of the joints of the system, the specification 
of displacement boundary conditions, the tabulation of 
material and section properties, the mass and stiffness 
matrices ace formed for each structural member and placed in 
sequence on low speed storage along with the global equation 
numbers wiich are associated with their stiffness terms. [In 
sa@aation, the member force-displacement transformation 
Matrices are formed and stored 9n a separate low speed 
storage file along with the appropriate displacement 
numbers. 

The NIDES operation is used td specify or generate the 
geometry of the system. The operation BOUND specifies which 
joint displacements exists and assigns internal equation 
numbers to these displacements. fTfherefore, each joint may 
have fron zéro to six displacezent degrees of freedon. 
Tables of material and section properties for the various 
menbers are loaded and print2ei as standard arrays of 
information. 

A sp2cial operation, : ADDSF, is used for the direct 
addition of element stiffnesses to form the total stiffness 
and diagonal mass matrix of the system. The ADDK operation 


may be used to add individual elements into the total systen 
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matrices. The LOADS operation specifies the concentrated 
joint loads for all load conditions. After the direct 
solution for joint displacements due to static or dynamic 
loads, the member fortes can be evaiuated useng the PORCE 
operation. Individual member forces can be evaluated using 
the MEMFRIO operation. The DISPL operation is used to print 


the displicements in joint number order. 





+ 
NODES ,41,41 
The cards following this Syogeaaeeee provide information for 
the creation of a arcay which will contain the 
coordinates for all the joi at's in the system. Where 

M1 


N1 = Number of joints (or nodes) in the systen. 


Name of new coordinate array to be loaded. 


The fh tale sequence of cards punched in a (215,3F10.0) 
format must follow this operation. 


Coluons Contain 
| ar. Node number selested by user 
6 - 10 blank 

11 = 20 X-coordinate 

24 = 30 Y-coordinate . 

31 = 40 Z~coordinate 


Nodes cards may. be supplied in IE order; however, all nodes 

most De dsfine2. Tf nodes are efined more than once the 
las ‘3 definintion will be used, fhis sequence of data must 
be terminated with a blank card. : 


Zz 
POSITIVE COORDINATE SYSTEMS 
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+ 
BOUND, M1 
This oOperition specifies the displacements which are nonzero 
for the structural system of joints specified by the NODES 
operation. Where 


M1 = Name of boundary condition cade array to be 
generated. 


This Qe OU ee followed by a_s2ries of cards containing 


the following information punched in a (815) format. 
Coluans Contents 
ie 5 Node number for the first node ina 


series of nodes with identical displace- 
ment specification. 


Ss = ie Node number for the last node in the 
series. : 

div= 25 X-translation. 

SG Y-translation 

Zee S Z-translation 

Boe 310 ) Co ao) oe toe ods! 

31 = 35 Y-rotation 

36 - 40 Z-rotation 

471 - 45 Node number increment used to generate 


Conditions for additional nodes. 


A translition or rotation equals: (a) zero for zero or 
undifined displacements, Or b) one for nonzero 
displacements to be evaluated by other operations. 


If a node boundary condition 1s not specified, all 
displacem2nts at that node are assumed zero. Cards may be 
Supplied in ay order. If node boundary conditions are 
Specified aore than once, the last difinition is used. This 
sejyuence of data must be terminatei by a blank card. 


The selection by th2 user of which nodes have nonzero 
displacements reguires an understanding of the direct 
stiffness procedure. Displacements degrees of freedom which 
have no stiffness associated with the displacement must be 
considered to be undifined since it 1S not possible to 
develop an equilibrium equation for that direction. The 
total number of nonzero displacements specified will be the 
Size of the total stiffness matrix to be defined by the 
ADDSF operation. 


116 





+ 
BEAM, M41,482,43,M4 


This operation calculates the el2ment stiffness, _mass and 
force-displacement transformation matrices for 3-D beam 
members. [These arrays are stored in sequence on low speed 
storage t>) be used by other operations where: 


M1 is the name of the beam element group 

M2 is the name of the coordinate array 

43 15 the name of the boundary condition array 

M4 15 the name of the arra which contain bean 
properties and has been loaded by the standard 
matrix LOAD op2ration 


On2 card for each beam in this group of beam elements nust 
follow this operation. The beam cards are punched in (515) 
format, where 


Coluans Contaila a. 
L oe Beam identification number 
Gr= 70 Node nunmber‘I 
Ue as Node number J. 
Ge 2.0 Node oumber K 
eae) 2 Beam property number NP 
This sequence of cards must be terainated with a blank card. 
The material and geometric properties for each element are 
give in the M4&Y array in the following order: 
M4 (N?,17) = Axial area of memb2r, A. 
MY (NP,2) = Torsional Moment of fnertia, J 
M4 (NP,3) = Moment of Inertia about axis 2, I 
MU (N2?,4) = Moment of Inertia about axis 3, I 
M4(N2,5) = Modulus of Elasticity, E 
M4 (NP,6) = Shear Modulus, G 
MU (NP,7) = Mass per unit length of bean 
where NP is t 


s the specific. material property number specified 
in columis 21 - 25 of the bean card. The local esacn 
convention 1S given in the following figure. 


aE? 





local beam reference system. 
Axes 1 and 2 are in the 
Plane defined by plane 

i, j and k nodes. Axis 1 

is defined by line i-j. 

Axis 3 is perpendicular 

the i, j, Kk plane. 





LOCAL SIGN CONVENTION 





DEFINITION OF POSITIVE BEAM FORCES 
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+ 
Pavia gl t,i2,03, 04 


This operation forms the element stiffness, _mass_ and 
force-displacement transformation matrices for 3-D truss 
- pembers. The arrays are stored on low speed storage in 
seyuence and will be used by other structural operations. 


M1 is the name of this groups of truss members 

M2 15 tne name of the coordinate array 

M3 15 the name of the boundary, condition array 

M4 15 an NP by 3 array of sestion properties in 
which NP is the number of different section 
properties and 


MU (N2?,1) = the cross-section area, A 
MG «NP,2) = the Modulus of Elasticity, 
an NP,3) = the mass per unit length of the member 


ration is followed by on2 card per truss member in 


rnoat with the following information: 

Coluaas Content ‘ cat . 

lL = .3 Truss member identification number 

6 - 10 Joint number I 

210 ody Res, Joint number J 

16 = 20 Section property number, NP 
This Operation must be terminated by a blank card. 
PS EAE gaa Po ri a 
LOADI,41,N1,N2,N3,N4 or LOADI,M1,N1,N2,,N4 Or 
EOADI,N1,N1,N2,N3 Or LOADI,41,N1,N2 
[This operation will load an integer ease named M1 which has 
N1 rows aid N2 columns. The terms of the array are punched 
in row-wise seguence on data cards which follow this 
operation. If 1S z2ro_ or. blank, the data must be punched 


in (1615) format. If N3 1s nonzero, an additiona card 
containing the format of the data cards must follow this 
Operation and receed the data cards. N4 is optional but 
note that an extra comma must be used in place of N3 if the 
user does not Supe a format. The matrix M1 will be 
pe aeee 1n partitione form with N4& columns per partition. 
ines hav> (N4 + 1)*5 characters. N& devaults to 20 causin 

125 characters to print per line. The user 1s cautioned no 
to overcome the capacity of the printing device in use. 
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. 
PLANE,M1,12,43,84,N1,N2 


This operation calculates the element stiffness mass and 
stcess-disolacement transformation matrices for 3 to 8 node 
lsoparametric elements. (Y-Z plane only). These arrays are 
stored in seguence aS a group on low Beeeo Betas to be 
used later by other operations (1.2., ADDSF and FORCE). The 
arguments are defined as 

M1 15 the user defined name of the element group 
M2 15 the name of the ae Ls coordinate array 
M3 is the name of the boundary condition arfay 
M4 is the name of the array which contains the 
material properties of the elements (one row per 


different material) where 


Not Modulus of Elasticity, E 
NP, 2 Poissons Ratio 
NP, 3 


Vv 
Thickness of the element 
= Mass density of the element 

NP is the material identification number 


Ni and N2 are the number of integration points in the r 
and s dir2actions respectively. 


On2 card for each 3 to 8 node element in the roup. must 
HoonoW ths Operation card. fTh= card is punched im a 
(1)I5,6F5.0) format and contains the following information: 

Columns Contain. rae 

lit 5 element identification number 

6 ~ 10 Node number N1 

11 - 15 Node numper N2 

16 = 20 Node number N3 

Pot) Ta I, Node number NY 

IAS © Bi Nodé number N5 

Slee Node number N6 

SiGe Node number N7/ 

Q1 - 45 Node oumber N8_. ; 

46 = 50 Material identification number, NP 

lin =O Natural coordinate? of stress output r1 

Oe 0 Natural coordinat?2 of stress output s1 

Stl 8 Tee) Natural coordinate of stress output 1r2 

66 - 70 Natural coordinate of stress output s2 

(2 is Natural coordinats of stress output r3 

76 = 80 Natural coordinate of stress output s3 
N4 throug1i N8 are optional. The nidside nodes, if present, 
must be #¢ithin center half of side. The local numbering 


System for the element 1S shown in the following figure. 








ISOPARAMETRIC ELEMENT 
Stresses will be printed by the PORCE a Aa at the three 


Olnts defined in columns 
efined as follows: 
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+ 
SLOPE, 41 


This operation forms a 4 x 4 stiffness matrix, M1 for_a bean 
or column member from the classical siope-deflection 
eguations. The properties of the member are defined on one 
card immediately following the operation. This second card 
1s punched in a (3F10.0) format and contains the following 
Infoegation: 


Columns Contain ; 

1=- 10 Moment of Inertia, I 
Les 220) Modulus of Elasticity, E 
21 - 30 Length of Member, 


The sign convention is defined as follows: 





DISPLACEMENTS 


| aS 


Mi 





V2 





M2 


L Mo | Ye 
FORCES FORCES DISPLACEMENTS 
The mMemoser forces are defined in terms of joint 


displacements by the following slope deflection equations. 


EI 
M,= — 





49, + 20, — a = “ 


E/ 
M, = Ee + 46, = $f» 


M,+M, 


Y,=-V = L 





Gr ihn watrix forn 


6 6 

4 2 L L 
M, ae 6 lie 
PA TE OW. lig. 
Vie i SG 2 2 lily 
V, | EO Oe be: L2{|r2 

& &§ _12 12 

L L L? 2 


or symbolically F = KA 
Matrix defined by th2 na 


a Ris) the 4 4x 4 stiffness 


+ + 
FRAME, 14,he 


This operation forms the 6 x 6 stiffness matrix 41 for the 
two-dimensional frame member shown below. 
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The properties of the member are defined on one card 
immediately following the FRAME operation card. This second 
card is punched in 3a (8F10.0) format and contains the 
following information; 


Coluans Contain 
1 - 10 Axial area, A a. 
11 - 20 Modulus of Elasticity, £ 
Zi = 30 Moment of Inertia, I 
31 - 40 Xi 
Gi - 50 vio 
51 - 60 Xj 
On ae O yes 
Up. 
Y Us; 
Uy; 
yj 
4) 


7 GEOMETRY AND 


7 
7 WA JOINT DISPLACEMENTS 
Aju 
a xt 


U 
y\ ; . — 


“1 x] < 


M2 is a3 x 6 force-displacement transformation matrix which 
1s based on the positive difinition of the element forces 


Shown below. 
a 
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These forses can be calculated from the following matrix 
equation, with the MEMFRC operation. 


U;, 
U,, 
M, Uy, 
M,; aa Af2 U,, 
Be U,; 
U 


+ 
LOADS, M1,482,N1 


This operation forms a load matrix named M1 of N1 columns 
(N17 load conditions) where M2 is the name of the boundary 
condition array Sha by .the operation BOUND. This 
Operation is followed by a series of cards - one for each 
loaded joint for each load condition. These cards are 
punched in (215,6F10.0) format as follows: 


Colusa 


S Contain 

5 Joint number 

0) Load condition nunber 
0 Load in X-direction 

0 Load in Y-direction 

0 Load in Z-directidsn 

Q Moment about X-axis 

: Moment about Y-axis 


n 
1 
Pu 
3 
4 
5 
6 
7 Soment about Z-axl1s 


OnUN & GIN) 


This series of cards must be terminated by a blank card. 
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+ + +¢ 
AUDSF ont on *'ADDSF,41, AZ 


[This operation forms the total stiffness matrix named M1 and 
a lumped aass matrix named M2 for the Structural system fron 
the element stiffness and mass matrices which are stored on 
low speed storade. These matrices can be printed with the 
PRINT opecation. If M2 is not specified, the row amass 
matrix M2 will not be formed. 


ADDK, 41,M2,83,N1 


This operation adds the element stiffness matrix named M2 to 
the total stiffness matrix named M1, where M1 was previously 
defined ani initially set to zero. M3 1s the name of the 
integer array in which the coluan number Nl contains the row 
Or column numbers in the total stiffness matrix where the 
element stiffness terms are to be added. 


L295 


o_o Se ee eins SO OY wR wp fe mee ee ee = es « °. c= eos; - —_ . 


2) gtd eld ".. average Ne wr eT A RTE, Cage ene SN eG tr Oe Oe ret nn Eee Tey! EP ane eee ee, CO3 oP MTS SL, Sars = 


oon we 





+ 
BESFRC,N1,82, M3, 44,N1 


This operation multiplies the element stiffness matrix named 
M1 by the joint displacement matrix named M2. 443 is the 
hare of the integer array in which the column number N1 
contains the row numbers, in the displacement matrix, M2, 
which are to be multiplied oy the element stiffness (or 
FoECe-disolacementy fhatrix, 1. Thegresuits Of HLS 
multplication are stored in the array named M4. 


This operation prints the. displasement named M1 in joint 
sejguence order, where M2 is the name of the bpoundary 
Comat tion array. 
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+ 
POmeCE Si, 82,43 or FORCE,N1;82 


This operation calculates the'member forces for a group of 
elements in whcih 
is the name of the element group 

M2 is the displacement matrix . ; 

M3 15 the name of the matrix in which the forces are 
stored ain the order calculated. GieniS A€AbrnayY 1S “Hot 
Specified, the element forces will be printed only and will 
not be retained in storage. For the TRUSS element, only the 
menber axial force, F, will be calculated for each member. 
For the 3EAM element, eight forces will be printed with 
reference to the positive definiton shown in the BEAM 
operation. 
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C. DYNAMIC ANALYSIS OPERATIONS 


The following operations were designed to evaluate the 
dynamic response of structures subjected to ap Ditrary 
tine-depenient loads. If these operations are used in 
connection with the standard matrix operations and the 
sturctural analysis operation, 32 dynamic analysis is a 
relatively simple procedure. The user has the option of 
uSing the sode superposition method or a direct step-by-step 
integration of the dynamic equations of notion. The user 
May examine the spectra of both input loading and calculated 
displacements. ine “aidicron, the contributions of the 
individual modes may be evaluated and compared. 

The 29st common and convenient form for time-dependent 
data to be specified is as straight line segments between 
given tim? points. Therefore, an operation which generates 
values at equal intervals is necessary. Another common 
characteristic of time-varying loads on structures 1s that 
it is normally possible to represent the loads at all points 
on the structure by the product of two matrices, a column 
matrix indicating the spacial distribution of loads times a 
Cow matrix which indicates the values as a function of 
various times. If a more complicated loading is required, 
it is possible to perform more analyses, each within the 
restrictions of the program, and then addthe results of 
each analysis. 

The following operations have been added for the major 


pucpose of performing dynamic analysis. 
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+ 
FUNG,M1,M2,M3,N1, N2 
This Operation generates a matrix named 82 which contains 
values at egual ‘intervals of the function specified in the 
eee ay named M1. The acray M1 must be a 2 by k array of the 
acm 


Cc 
2 os, 


ae numerically represents a function of the form shown 
below: 


M1 = 





tf, 


The tine interval t is specified in the 1x 1 matrix tee 
N1 specifies the Oran number of values to 
ipecarrsy" and is the number of columns in M2. N 
he arra un wee be a 1x Nl row matrix in which the a 
value wi fe I= N2 O, the array M2 will be a 2 x N1 
matrix a tbs following obi: 


ty t, + At t + 2At #eeeeeee##e#e#8e#e#e2s @ 





f f(t, +At) f(t, 4 Qt) oes 8 a 
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- - + 
Soiree), M2, ha, MS Moo Ay , AS,NI,N2 


This operation calculates the dynamic response (of a 
sturctural systen Se direct Se P= Eiete P Lowegration of 
-~the following linear matrix equation of motion: 


MO + CU + KU = R(t) = PF(t) 


Where 
M1 is the name of the N x N stiffness matrix K 
M2 15 the name of the N x N mass matrix # 
M3 is the name of the N x N damping matrix C 
M4 is the name of the N x 3 initial condition matrix 


0 in which: 


0 (I,1) is a vector of displacements JU 
UM tSea Vector Of velocities J. 
U 1s a vector of accelerations JU 





ie 
M5 15 fhe name of the N x N2 gatrix of calculated 
displacements in which column "i" represents the 
displacements at time 1i*N1*At ; 
M6 - he name of the N x:'1 load distribution matrix 


M7 15 the name of the 1 x k row_matrix representing 
the load anh tr ial at equal time increments 
FP, where k = N2/N1 : ; ——- 

M8 is the name of the 1 x 1 matrix containing at 

N1 15 the output interval for the displacements 

N2 is the total number of displacement vectors to 
b2 calculated. 


The total time for which results will be calculated by this 
Op2ration 15 Ni*N2*At. This operation must be followed with 
one data card in (3F10.0) format containing the following 
1nformation: 


Colunas contain 
1 - 19 DELIA 
11 - 20 ALPHA 

21 - 30 THETA 


Different values of delta, alpha and theta will allow the 
USer to select different methods of Shee -step 
lntegratijn. The following table lists some possibliities: 


DELTA ALPHA THETA 


c 


Newmacks Average Acceleration V2 1/4 0 
Linear Asceleration 1f2 1/6 1.0 
Wilson's Theta Metnod (low damping) 1/2 1/6 1.42 


Wilson's Theta Method {high damping) 1/2 1/6 2.0 





- + = 
ProiN,M i,h2,n3, Nt 


This operation solves the following eigenvalue problen: 
Kg = MGA 


In which the Nx N, symmetric, positive-semidefinite matrix 
K is named M1. the matrix M 1s a diagonal matrix of 
nonzero, positive. terms designatei by M3. Tha @atrix 3 
must be a row or column matrix containing only the diagonal 
terms of 4M. The e€igenvalues,A, are stored in matrix M3. 
Th? eigenvalues are ordered in hunerically increaSing order 
and the 2igenvectors,g@, are stored in the corresponding 
columns o£ the matrix M2. The number Nl specifies the 
approximaze number of Significant figures of the 
2ijgenvalu2s. If Ni is zero or blank, 4 figure accuracy will 
be used. fhe maximum accuracy possible is 16 figures. The 
us2 of mor2 than 12 figure accuracy is not recommended. 


The ee ce reduces th2 problem to standard eigenvalue forn 
+ 


by t Ollowing transformation 
where 

Ky = mt Km 

I =a Mo 
mm which 

m; = L/ /May 


The calculated mode shapes,@, are normalized as follows: 
f'ué = 1 g*KG A 


The program uses the standard Jacodi diagonalization method 
to solve for all eigen values and 2igenvectors. 
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+ 
DYNAM,M1,42,M3,M4,M5, M6 ,N1 


This operation evaluates the following set of uncouples 
second order differential equations associated with the node 
Superposition method for the dynamic analysis Om Ga 
Structural systen. 


ee e a ms : 

Xj + 2h; Wz X; el Aue Fae) i = 1 to N nodes 
M1 is the name of a row or column matrix which contains the 
N terms de a in rad/ses). M2 1s the name of a 
crow or column matrix which contain the N A; terms (ratio of 


modal damping to critical damping). 


The generalized time-varying forces P; (t) are not specified 
dicectly but are evaluated fron more fundamental 
information, The forses for all modes are _ evaluated at 
Specific times by th2 program from the following matrix 
equation: 


P = p*f = M3*M4 
x 


In which p 1S a specified N ‘1 vector named M3, and isa 
ee GOW Macelx which wiidl? be ere eaec from the by ok 
4 1s he same form as the input 


arccay nam2i M4. The Lee M4 i 
arcay described under he operation FUNS. rt 15 Hot 
neceSsary to use FUNS before the DYNAM operation. 


M5 ais the name of the N x Ni array which contain the 
generalized displacement X; (t). 


M6 is the name of the 1 x 1 array which contains the time 
increment associated with the generalized displacements. 


N1 1s the number of displacements to be generated. 


The method of integration used is exact for straight line 
secments. 


Wey) 








PLOT, M1,N1 


This operation will Suet ieic a printer plot of selective rows 
of the matrix named M1. N12 is the number of rows of M1 
wulchewill betplotted by ™Ehis operation. This operation is 
followed by Ni cards in (iat, raf format with the following 


information: 


Coluans Contain 
1 Plot symbol - any each symbol 
Leo Row number to be ploted 


The pr Odrag automatically searches the information to be 
plotted f5r the maximum and minimum values. The difference 
an aera numbers diviied by 120 spaces is selected as the 
pikgt SCai2. 
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D. LOOPING OPERATIONS 


CAL has a five level looping ability. The first 
Operation is LOOP and the last operation is  WEG?. 
Operations withing CAL are normally executed as they are 
encountered. If the operation reguires data, the data cards 
follow the operation card. In the case of looping, however, 
all operation cards from the first LOOP card to the last 
NEXT card are stored within the computer before they are 
executed; therefore, if operations within the loops regure 
data, the data cards must be, supplied in the order required 
after the last NEXT operation. If an error is encountered 
while ex2cuting in a loop, tha entire smatrix of loop 
commands is deleted and the user is given the opportunity to 
try again. Matrices that have been modified by operations 
successfully completed while in the loop remain modified. 
After all loops are executed the computer storage required 
for thes? operations iS automatically released by the 


program. The looping operations are: 


hoe 





LOOP, N1 


N1 ais the number of times the loop is to be executed. 
Associatei with each LOOP operation there must be a 
corresponding NEXT operation which signifies the end of the 
Loop and the return of the control,.to the beginning of the 
loop. . Tae following is a possible series of looping 
Operations. 


LOOP, 5 














First level loop 


LOOP, 2 


second level loop executed 5 times 


executed twice 
(total of 10 times) 


NEXT 


LOOP, 4 Second level <loop 
oe Jekeieve executed 4 times 
NEXT (total 20 times) 
NEXT: 


NEXT,M1 or NEXT 


The operation NEXT signifies the end of a loop. It is 
apparent which LOOP and NEXT caris are associated if there 
are an egual number of each. The operation _NEXT,S1 will 
cause the loop to terminate if the first term in the matrix 
nazned M115 negative. 


SKIP,M1,N1 

This operation will caase the skip of the next N1 operations 
if the first term in the matrix named M1 1S negative. This 
level of looping. 


eS me eee 2 Ss See 2 SP we @S @ ew SS ws ae se SP Sew a ee eS Se Se ew ee Se | Se Se Se ee ee Ss we 6S | 6 2 Se 6 So = GF = 
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E. USER DEFINED OPERATIONS 


USERA and USERB 


These nam2s are reserved for operations to be defined and 
programmei by the 1ser. In” order to program these 
Operations 1S 1S necessary tod understand the internal 
organization of CAL. Chapter III contains details. 
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¥. LARGE PROBLESS 


CAL 1s designed as an educational tool. It does not 
take advantage of banding and symmetry in matrix storage. 
Larger problems can be solved by increasing the dimension of 
the L array but a general purpose program that makes maxinun 
advantage of out-of-core storage and takes advantage of 
banding and symmetry for in-core agjatrix storage is probably 
a better choice. With the above disclaimer, to increase 
problem siz2 capability, increase the dimension of the L 
array and change the value of MAX to the new dimension size 
wee ch LOL losing: | 

Cae MAIN PROGRAM 
=== SET PROGRAM CAPACITY 
COMMON NTOIT,NDP,L(6000) 
MTOT=5 000 
ND P=2 
SAL. 5 ELEN E 
SALE CALI 
STOP 
END 
SETIME is an OS routine to initialize the CPU tiner and 
Should not be used for CP/CHMS programs. With the dimension 
of the L array as above, the program currently executes in 
144k bytes for OS and 256k bytes for CP/CMS. The region 
necessary for execution will increase about eight times the 


increase in the L array. 
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